On the conservation of spin currents in spin-orbit coupled systems 



in 
o 
o 

(N 

o 

Q 

in 



R. Shen, 1 - 2 Yan Chen, 1 and Z. D. Wang 1 - 2 

1 Department of Physics, The University of Hong Kong, Pokfulam Road, Hong Kong, China 
2 National Laboratory of Solid State Microstructures and Department of Physics, Nanjing University, Nanjing, 210093 China 

Applying the Gordon-decomposition-like technique, the convective spin current (CSC) is extracted 
from the total angular-momentum current. The CSC describes the transport properties of the 
electron spin and is conserved in the relativistic quantum mechanics approach where the spin- 
orbit coupling has been intrinsically taken into account. Arrestingly, in the presence of external 
electromagnetic field, the component of the convective spin along the field remain still conserved. 
This conserved CSC is also derived for the first time in the nonrelativistic limit using the Foldy- 
Wouthuysen transformation. 
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Spintronics is a fast developing field in condensed 
matter physics, where the research on the spin transport 
and manipulation in solid state systems appears to 
be a central subject ^£ Typically, the spin transport 
is strongly affected by a coupling of spin and orbital 
degrees of freedom in semiconductors. Of particular 
interest is the spin Hall effect in a spin-orbit coupled 
system. The spin Hall effect means that a trans- 
verse spin current can be induced by applying an 
electric field. Earlier studies focused on the extrinsic 
effects due to the spin dependent scatteringi 3 i 4 i 5 Re- 
cently, an intrinsic spin Hall effect (1SHE), which is 
resulted from the spin-orbit coupled bands and may 
exist in the clean systems, has attracted considerable 
attention: 6 - 7 - 8 - 9 - 10 - 11 - 12 - 13 - 14 - 15 - 16 - 17 - 18 - 19 - 20 - 21 - 22 - 23 - 24 - 25 - 26 - 27 
On the theoretical side, the 1SHE has been predicted in 
the p-doped semiconductors described by the Luttinger- 
like models and in the two-dimensional electron gas 
described by the Rashba-like model 7 . The 1SHE in the 
two-dimensional hole gas^Si and the impurity scattering 
effect on the ISHE 15 i 16 i 17 i 18 i 19 i 20 i 21 have also been ana- 
lyzed. On the experimental side, the spin accumulation 
related to the spin Hall effect has been observed in the 
n-type bulk semiconductor— and in the two-dimensional 
hole gas— 

The spin current is an important issue in discussing 
the spin Hall effect and other spin transport phenomena. 
In spite of the extensive investigation on the ISHE, there 
is still a fundamental problem that what is the proper 
definition of the spin current operator in the spin-orbit 
coupled system. The conventional definition of the spin 
current operator is the anticommutator of the Pauli spin 
matrix and the velocity of the electron^ 7 . The electron 
spin is not a conserved quantity in the spin-orbit cou- 
pled system. As a result, there is a source term, namely 
spin torque, in the continuity-like equation related to the 
conventional spin current— The spin torque repre- 
sents the spin precession caused by the external field and 
the spin-orbit coupling. The presence of the nonzero 
source term in the continuity-like equation makes the 
conventional spin current an ill-defined quantity to de- 
scribe the transport of spin. On the other hand, the 
non-conserved part of the conventional spin current has 



very fast dynamics and can not contribute to the low fre- 
quency response— Therefore, it is essential to define a 
conserved effective spin current from a first principle for 
studying the spin transport properties in the spin-orbit 
coupled system. 

There are already some efforts focusing on the under- 
standing of the nonconservation of the conventional spin 
current i 22 i 23 i 24 i 25 i 26 i 27 i 28 i 29 In the p-doped semiconductor 
described by the Luttinger-like model, Murakami et al— 
have proposed a conserved effective spin current of the 
spin 3/2 hole by modifying both the spin density and 
the spin current operator. By rewriting the spin torque 
as the divergence of a torque dipole density, another ap- 
proach to the conserved effective spin current was also 
28 ^^SiedSML where the spin density operator remains 



unchanged and only the spin current operator is modi- 
fied. In this paper, motivated by the fact that the spin- 
orbit coupling in the semiconductors comes from the rel- 
ativistic correction of the electron motion, we develop a 
conserved effective spin current in the relativistic quan- 
tum mechanics approach where the spin-orbit coupling 
has been intrinsically taken into account. For a rela- 
tivistic electron, the spin is not conserved, but the total 
angular-momentum (TAM) consisting of both spin and 
orbit angular-momentum is a conserved quantity in the 
absence of external fields. Therefore, there is the con- 
served TAM current according to the Noether's theorem. 
Under the Gordon-decomposition-like scheme, the TAM 
current can be divided into two parts: the convective 
TAM current and the internal one. It is found that the 
spin part of the convective TAM current, namely the con- 
vective spin current (CSC) satisfies the continuity equa- 
tion. The CSC is not only conserved but also closely 
related to the spin flux of the electron moving in the spa- 
tial space and therefore is an appropriate candidate for 
the spin current operator in discussing the spin trans- 
port properties. In the presence of external electromag- 
netic field, neither the TAM nor the CSC is conserved 
in a general case because of the mutual exchange of the 
angular-momentum between the electron and the exter- 
nal field. However, more intriguingly, we find that if both 
the electric field and the magnetic field are along the same 
direction, namely the z-direction, the z-components of 
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the CSC and the TAM still remain conserved. Remark- 
ably, this is just the case in the Rashba-likc spin-orbit 
coupled system. As a result, the CSC is also likely ro- 
bust subject to the external field in some sense. Finally, 
by using the Foldy-Wouthuysen transformation, the con- 
served CSC is derived in the nonrelativistic (NR) limit, 
which has a potential application in spintronics. 

We start with the brief review on the Gordon decom- 
position of the probability current of the relativistic elec- 
tron. For the electron of the mass m traveling in the 
three-dimensional free space, the four-component wave 
function of the electron m satisfies the Dirac equation 



dm 

at 



(cot • p + mc 2 0) m , (1) 



where the momentum p 



(8 



dfi,V, the matrices a 



a is the Pauli spin matrix and 1 



is the 2x2 unity matrix. The invariance of the Dirac 
Eq. under the U(l) gauge transformation results in 
the continuity equation for the probability current 



(2) 



where the probability density p — ty'ty and the current 
j = cm 1 ' am. The Dirac Eq. can be rewritten as 



mc 



1 ./.J h- 
—p [ in— — cot* — V 



dt 



i 



m. 



(3) 



Substituting Eq. © into the definitions of p and j, one 
finds that the density and current can be divided into 
the convective and internal parts, known as the Gordon 
decomposition^SiS in the form 



P = Pc 



Pi, 



3 = Jc + Ji, 



(4) 



rate of change of m in space and time and closely related 
to the motion of the electron in the spatial space. Notice 
that the convective and internal parts of the density are 
separately conserved. For the NR limit, p x is zero up to 
the zeroth order in (1/c), j\ is only a curl term and has 
nothing to do with the electron accumulation, and the 
convective parts of the current and the density are re- 
duced to the usual form in the NR quantum mechanics, 



(V^)V], (8) 



with ip the two-component wave function of the electron. 

Next, we turn to the spin current. According to the 
Noether's theorem, the invariance of the Dirac Eq. JJJ 
under the Lorentz rotation leads to the continuity equa- 
tion for the TAM current 

^ + V.j 7 = 0, (9) 

where the density of the 7-component of the TAM is 



P 1 = \{^J^ + h.c.}, 

and the current of the 7-component of the TAM is 
1 



(10) 



ft = -c{*W 7 * + /»•&}, 



(11) 



where the 7 component of the TAM operator = 
((7i/2)S + r x p) 1 with the position operator r and the 
index 7 = x, y, z. The TAM contains two parts, the spin 
part (7i/2)E and the orbit part r x p. For a relativistic 
electron, unlike the TAM, neither the spin density nor 
the orbit angular-momentum density is conserved alone. 
Similar to the probability current, the TAM current can 
also be divided into the convective and internal parts^ 



2mc 2 



m 



om dm 

~dt~~dt 



m 



2mi 



(5) 



2mc 2 



dm dm 

~dt~~dt 



h.c. 



2 I2^i [*V(J 7 *)-(V*)J 7 *] 



(12) 



V • P, ji = cV x M 



dP 



where 



(6) 



m = — mum, p = —m{-ia)m, (7) 

2mc 2mc 

and the matrix £ = ( £ ° ) as wen as W = m^fi. The in- 
ternal parts of the density and current are closely related 
to the internal state of the electron. For the electron of 
charge q, qM is essentially the magnetic moment density 
of the electron arising from its spin and qP is the elec- 
tric polarization belonging to the spin. The convective 
parts of the density and current are determined by the 



Pi 

where 



V«F 7 , j 7 = cVxM 7 + 



Qp~t 

w 



(13) 



M =l( J-$EJ 7 f + h.c. I , 



2 \2mc 



37 



1 



2mc 



m (-ia) + h.c. 



(14) 



The convective and internal parts of the TAM current 
have the same structure as those of the probability cur- 
rent. For the NR limit, p 7 is zero up to the zeroth order in 
(1/c), f is only a curl term and has nothing to do with 
the accumulation of the TAM. Let's consider the spin 
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part of the convective TAM current, namely the CSC, 
with the density for the spin-7 component 



n s "f 



xh 2 - 

4mc 2 7 dt 



h.c. 



and the current 



~s-7 _ 



2mi 



*^£ 7 (V*)-(V*)^ 7 * 



(15) 



(16) 



The CSC is determined by the rate of change of ^ in 
space and time and ought to be closely related to the spin 
flux of the electron moving in the spatial space. More- 
over, noting that the solution of the Dirac Eq. JJ) is also 
a solution of the Klein-Gordon equation, one can easily 
find that the CSC satisfies the continuity equation 



dp, 



,S"7 



dt 



v.jT 7 = o, 



(17) 



despite that the bare spin is not a conserved quantity 
for the relativistic electron. Therefore, we argue that 
the CSC is just the correct effective spin current for dis- 
cussing the spin transport properties of the relativistic 
electron. 

The results we present above are valid for the electrons 
moving in the free space. The influence of the external 
electromagnetic field can be taken into account by mak- 
ing the transformation 



P 



c dt 



d 



(18) 



with the charge of the electron q, the scalar potential V 
and the vector potential A. The density and current of 
the convective spin in the external electromagnetic field 
are written as 



/ ih 2 — d^ 



■>■-"■■ - I r#£-v — + h.c. ) - -qV^Z-y^, (19) 



\4mc 2 



dt 



2mc 2 



•s-7 



— ^ l{ V*)- h ,]-«—^. (20) 



In the external electromagnetic field, neither the TAM 
nor the CSC is conserved generally because of the mutual 
exchange of the angular-momentum between the electron 
and the external field. The continuity-like equation of the 
CSC reads 



dp: 



,s-7 



dt 



V.ir 7 = /* B *(SxB-iax£) % (21) 



where the Bohr Magneton /ib = qh/2mc, the mag- 
netic field B = V x A, and the electric field E — 
-(dA/dt) - W. Although the CSC is not conserved, 
the source and drain of the CSC in Eq. (|21|l are only the 
external fields, while those of the bare spin include the 
contributions from both the external fields and the or- 
bit angular-momentum. Therefore, the CSC is still more 



suitable than the bare spin current for discussing the spin 
transport of the relativistic electron in the external fields. 
For a particular case that both the electric and magnetic 
fields are along the same direction, namely the z direc- 
tion, one can find from Eq. (|21() that the z-component 
of the convective spin is still conserved. This particular 
case is relevant to the two-dimensional electron gas with 
the Rashba interaction where the electron is free in the 
x-y plane and the confinement potential is only along the 
z direction. The z-component of the spin is also of the 
most interest in the Rashba system. In the following text 
we are restricted to this particular case, where the CSC 
is still conserved and can be served as the proper quan- 
tity for discussing the spin transport of the relativistic 
electron. 

In order to apply the idea of the CSC into the spin- 
tronics, the NR limit of the CSC operator is needed. The 
Hamiltonian of the relativistic electron in the external 
electromagnetic field reads 



H = cot 



+ qV + mc 2 /3. 



(22) 



The NR limit of the Hamiltonian (1221) can be made by the 
Foldy-Wouthuysen transformationi^iSi Up to the (1/c 2 ) 
order, the NR Hamiltonian of the electron is given as 



1 qti 
H NR =— 7T 2 + qV ~ p B cr • B - 

Zm &m z c z 



V *E 



qti 



8m 2 c 2 



er • (E x n — ir x E) 



1 



(23) 



j Pb 

8m 3 c 2 2mc 2 

Mb 
4m 2 c 2 



(E 2 - B 2 ) 



[(<T • B) 7T 2 + TV 2 (CT.B)] , 



and the two-component wave function ip, i.e., the 
large component of the four-component spinor ^ in the 
Foldy-Wouthuysen representation corresponding the pos- 
itive energy solutions, satisfies the Schrodinger equation 
ih(dtp/dt) — iJNR/0- The Schrodinger equation and the 
Hamiltonian l|23|) are the starting point for discussing the 
motion of the election in the NR quantum mechanics in- 
cluding some relativistic corrections up to the (1/c 2 ) or- 
der. We recall the physical meanings of each term in Eq. 
(|23|l . The first three terms are the NR kinetic energy, the 
potential energy, and the Zeeman energy arising from the 
electron spin, respectively. The fourth and fifth terms are 
the so-called Darwin term and the spin-orbit interaction, 
respectively. All other terms are of even higher orders. 

With the help of the Hamiltonian 1221) and the Dirac 
equation, the convective spin density in Eq. (|19|) can be 
rewritten as p s ~ z = p s c ~ Zl $> + h.c.) /2, with the opera- 
tor 



2 2mc 



(3 {ax 7r) 2 



(24) 



One can easily verify that p s c ~ z is an Hermitian opera- 
tor and commutes with the Hamiltonian 1221 and there- 
fore is a conserved quantity. The difference between the 
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convective spin and the bare spin is the second term 
in Eq. (|24[l which is a higher order relativistic correc- 
tion term coupling the large and small components of 
the four-component spinor 'J'. The NR limit of the con- 
vective spin density operator up to the (1/c 2 ) order, the 
same order as that of the Hamiltonian, is obtained by the 
same unitary transformations applied to the Hamiltonian 
(|22|1 in the Foldy-Wouthuysen scheme. After the Foldy- 
Wouthuysen transformation, the convective spin density 
is also reduced to an even operator as the Hamiltonian 
does, which does not couple the large and small com- 
ponents of \P. The NR convective spin density operator 
acting on the two-component spinor ip is obtained as 



with the current 



Pc,nr — 



2 



4m 2 c 2 

n 

Am 2 c 



(<T • 7f) (<T X 7r) 



-c. 



(25) 



with the operator C = (it 2 . + 7r 2 ) a z — (<r x ir x + <Tyir y ) tt z . 
With the help of the operator l|25|) the convective spin 
density is written as P^nr = (^PcNrY 1 + ^- c -) /2- Due 
to the unitarity of the transformation, the NR convective 
spin density operator Eq. (|25() commutes with the Hamil- 
tonian and remains as the conserved quantity up to 
the (1/c 2 ) order. The difference between the convective 
spin and the bare spin is the second term in Eq. 125|) 
which is of the order of (1/c 2 ) and is very small. How- 
ever, the spin-orbit interaction is also originated from 
the relativistic correction of the electron motion up to 
the (1/c 2 ) order. Therefore, the second term in Eq. I|25|) 
can not be neglected when discussing the spin current in 
the spin-orbit coupled system. 

Although there are many terms in the NR Hamiltonian 
(|2lj) . only the spin-orbit term is of particular interest on 
the topic of the spin Hall effect. Therefore, we remove 
the other relativistic corrections and rewrite the effective 
Hamiltonian in the form 



Hso 



2m 



7T + qV - [IbB(t z 



qhE 

Am 2 c 2 



(<TX7r) z . (26) 



It should also be noted that, in going from Eq. (12.'j[) to Eq. 
(|26|l . it has been assumed that the electric field is static, 
i.e., E = -VV, and the external fields are along the z- 
direction. It can be easily verified that p^^ R commutes 
with -ffso and therefore served as a conserved quantity to 
calculate the spin flux in the spin-orbit coupled system. 
With the help of Eqs. 1 (25) 1. J2EJj and the Schrodinger 
equation \h(dtp/dt) = Hsoip, the continuity equation for 
the CSC in the NR limit is obtained in the form 



dp, 



c.NR 



dt 



(27) 



Jc.NR 



4to 



h.c 



(28) 

Eqs. Q25[). (|27|l . and l|28|l constitute the main results of 
this paper. Comparing to the conventional spin currentji 
the CSC (|28|l has an additional correction term 



16m 3 c 2 



(29) 



which comes from the second term in Eq. (|25|) . One 
finds that both the spin density and current have been 
modified to construct the conserved CSC. The modifi- 
cation, the second term in Eq. I|25|l and Eq. I|29|) . is 
a kind of the relativistic correction resulting from the 
coupling between the large and small components of the 
four-component spinor <J . Those corrections are of the 
(1/c 2 ) order and are often small enough to be neglected. 
However, the spin-orbit interaction is also a relativistic 
correction of the (1/c 2 ) order. Therefore, the CSC ought 
to be used instead of the bare spin current in the spin- 
orbit coupled system. 

Although the Hamiltonian (|26|) looks similar to the 
Rashba model2L2& in the semiconductor heterostructure, 
a realistic situation is quite complicated. The Rashba 
coefficient a may not be directly determined by the ex- 
ternal electric field, but the asymmetry in the band struc- 
ture parameters of the heterostructure like the effective 
massi 2 *^ The discussion in this paper is restricted to the 
case that the electron moves in a uniform space. There- 
fore, it is interesting and desirable to apply the above idea 
of the conserved CSC to the semiconductor heterostruc- 
tures or quantum-wires in future, though is highly non- 
trivial. 

In summary, the CSC has been extracted from the 
TAM current by the Gordon-like decomposition in the 
relativistic quantum mechanics approach. The NR limit 
of the CSC is also derived. The CSC describes the trans- 
port properties of the electron spin and is conserved in 
the spin-orbit coupled system. It is also shown that the 
CSC remains robust in the external electromagnetic field 
if both the electric field and the magnetic field are along 
the same direction. As a result, the CSC is a good sub- 
stitution for the bare spin current in studying the spin 
transport phenomena in spin-orbit coupled systems. 
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